Non-adiabatic coupling and adiabatic population transfer 
in quantum molecular systems 
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We show that a counter-intuitive pulse sequence leads to adiabatic passage between the vibrational 
levels of three harmonic potentials through parallel dark states in adiabatic approximation. However, 
the adiabatic assumptions break down for very intense pulses and non-adiabatic couplings result in 
the population transfer by light-induced potential shaping. 
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Population transfer by shaping light induced potentials 
(LIP) is a very robust scheme to transfer vibrational pop- 
ulation between different electronic potentials in diatomic 
molecules 0, |[ ||, |[ ||, |j. One of the main advantages 
of the scheme is that the vibrational quantum number 
is conserved during the overall process. Since the selec- 
tivity of the transfer is guaranteed by a symmetry rule, 
the method is not very sensitive to the exact position of 
the energy resonance, i.e. to the topological or energetic 
features of the potentials involved. 

The simplest scenario of population transfer by shap- 
ing LIP can be described in terms of three electronic 
states (potential curves) sequentially coupled by two 
laser pulses, and the vibrational population is driven by 
two photon off resonant absorption. Although there are 
several possible scenarios for selective population trans- 
fer (HQ, the first scheme proposed, called APLIP, in- 
volves a sequence of two very strong (TW/cm 2 ) picosec- 
ond pulses applied in counter- intuitive order The 
pulse that couples the intermediate electronic state with 
the final excited electronic state must precede the pulse 
that couples the ground electronic state with the inter- 
mediate one. Then the system dynamics follows a char- 
acteristic pattern in the vibrational basis representation, 
which is illustrated in Fig. 1 for a test system of three 
harmonic oscillators with the same force constant. (More 
details about the model are presented below.) According 
to Fig. 1(a), the overall population on the ground elec- 
tronic state is rapidly transferred to the final electronic 
state, while very few population is temporally excited to 
the intermediate electronic state. So, in principal it is 
possible to suppose there is a dark state consisting of 
the initial and final potentials, similar to the population 
transfer through the dark state of the three-level system 
H . Looking into the vibrational populations in detail, it 
can be seen that the passage is mediated by substantial 
excitation of high energy vibrational levels in both initial 
and final potentials, Fig. 1(b). 

The scenario previously considered is very similar to 
the well known Stimulated Raman Adiabatic Passage 
(STIRAP) scheme g. In STIRAP there are three levels 
sequentially coupled by two fields working in counter- 



intuitive order. This arrangement of pulses prepares the 
system initially in a dark state, |$o)j which is an eigen- 
state of the Hamiltonian dressed by the field in the adi- 
abatic approximation. The dark state correlates at late 
times with the final level |3) and never overlaps with the 
intermediate level. The population is transferred from |1) 
to 1 3) in a completely similar way to the overall electronic 
populations of Fig. 1(a). 
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FIG. 1: Population dynamics in three symmetrically dis- 
placed harmonic oscillators model, (a) Solid line shows the 
total population of final state, dashed line shows the total 
population of ground state, and dotted line is the population 
of intermediate state, (b) Population of vibrational levels, 
dashed line - zero vibrational level of ground state, dot-dashed 
line - population of the vibrational levels v > of ground 
state, solid line - population of zero vibrational level of final 
state, dotted line - population of the vibrational levels v" > 
of the final state, (c) Shape of the laser pulses. 

Normally, the STIRAP scheme involves nanosecond (or 
longer) pulses with intensities in the GW/cm 2 (or less). 
Therefore it is natural to ask if APLIP is a short time, 
strong field version of STIRAP (and applied to electronic 
states instead of vibrational or atomic levels). Indeed the 
requirement of stronger laser fields can be understood un- 
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der the general assumptions of STIRAP. Since the adi- 
abaticity condition in STIRAP is usually expressed by 
the pulse area relation, fir 3> 1, where f2 is the effective 
Rabi frequency of the pulses and r their time widths, any 
reduction in the pulse width must imply a correspond- 
ing increase in the pulse amplitude. Nevertheless, in this 
note we show that the population transfer in APLIP can 
not be reduced to the generalized description of STIRAP 
in multi-level system. In order to do so we show that 
the natural extension of STIRAP to the three potential 
curves scenario in the strong field limit predicts a differ- 
ent dynamical evolution than the one shown in Fig. 1. 



FIG. 2: (a). Scheme of three symmetrically displaced har- 
monic oscillators truncated to give two coupled 3-level ladder 
systems, (b). Dressed potentials, Ui = Vi, U% = V% — wi, and 

Us — V3 — U)i — L02. 



Let us consider in detail the Hamiltonian of the system. 
The Hamiltonian for three electronic potentials (Vi(x), 
i = 1,2, 3) coupled by two laser fields (-Ei(t) and E 2 (t)) in 
the Born Oppenheimer and rotating wave approximation 
(RWA) reads 



T + U ± (x) -^ 12 {x)E 1 (t) 
H RWA = [ -^ 12 (x)E 1 (t) T + U 2 (x) -y 23 (x)E 2 (t) , (1) 
-±fi 23 (x)E 2 {t) T + U 3 (x) 



where T is the kinetic energy operator and Ui(x) are 
the diabatic dressed potentials (U\(x) = Vi(x), U 2 (x) = 
V 2 (x) — u>\, U 3 (x) = V 3 (x) — lui — uj 2 , LOi being the car- 
rier frequencies of the lasers), ^i 2 (x) and ^ 23 (x) are 
the dipole moments (atomic units are used through- 
out). In this representation the global wave func- 
tion of the system is expanded in terms of electronic 
wave functions, {B.i(q; x)}, and vibrational wave packets, 
{ipi(x,t)}, where q is an index representing the collective 
electron coordinates and x is the vibrational coordinate. 

After expanding the vibrational wave packet in the set 
of eigenf unctions of the vibrational Hamiltonian for each 
potential, 

i; a (x,t) = J24 a) W<t>T ) W> ( 2 ) 
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where (%) is the j'th vibrational level in the electronic 
state a (in the following we use Greek letters to designate 
electronic states and Roman letters to designate vibra- 
tional levels), we arrive at the following time dependent 
Schrodinger equation (TDSE): 

< id? = -T k f4 ) + {^ -Dl) d ? - Y, k ^df 

. 4 3) =-£* ?4 2) + {»? - -2 + Dl Dl)df , 

(3) 



where u>j + D? a -\ is the eigenvalue corresponding to the 
\(j>f ] ) eigenfunction, fig 5 = E 1 {t){4>f ) ^ 1 \ i i 1 \4>f ) ^ 2 ) and 
Q\j = E 2 (t)(4>\ E 3 \fi 2 \(j)j 'E 2 ) are the Rabi frequencies, 
and are the potential zero energies. 



In order to compare the dynamics of the general sys- 
tem with that of STIRAP, we establish a correspondence 
between the 3 x N equations of motion (Eq. (||)) and 
the equations for N 3-level ladder systems, where both 
diabatic and adiabatic states are known, corresponding 
to that of the STIRAP Hamiltonian. In doing this con- 
nection we neglect all the contributions from the con- 
tinuum wave functions in all the electronic states. To 
simplify the notation we use a model of three symmet- 
rically displaced harmonic oscillators (Fig. 2) with the 
same force constant (SDHO model), so that the sep- 
aration between the minima of the potentials is con- 
stant, 7q 3 ^ — r$ = 7q — = R. We assume that 
the equilibrium configurations of the excited potentials 
are displaced to larger inter nuclei distances. Further- 
more we consider only processes in two photon reso- 
nance. This allows us to define a constant energy split- 

4.- X (!) (!) (3) (3) j i i 

ting, duo — oJ^+i ~ w « — k-n-i-l ~ w ™ i an d a one photon 
detuning, A = L>§ - c^ 1} - wt = (oj 2 + D%) - uj[f ] - D%. 
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The TDSE is: 



il<l? ,(2) 



3 

j(3) _ 
l 3 ~ 



2 U k 



,(2) 



,(2) 



Wd 3 

-E,^4 2) + (^-A)< 



(3) 



(•'!) 



(4) 



Finally we invoke the Condon approximation. There- 
fore, the state to state Rabi frequencies can be expressed 
as fig> = pjfuxExit) and = pf?/i 2 S 2 (i), where 

we have defined the geometric Franck-Condon parame- 
ters p^ 3 = (<^>i |<^j )• This is the model used to ob- 
tain the results shown in Fig. 1, where we have cho- 
sen Sui = 7 • 10~ 4 a.u., A = 0.015 a.u., \i\ = [i 2 = 1, 
E l (t) = E S(t), with E = 0.08 a.u. and S(t) being a 
Gaussian envelope function. Both pulses have the same 
envelope function, with width a = 2.5 ps and they are 
time delayed, so that the second pulse precedes in 2.5 ps 
the first pulse. To obtain the numerical results shown in 
Fig. 1 we have numerically integrated the TDSE based on 
a grid discretization of the Hamiltonian in Eq. (|l]) and not 
on the discrete basis representation of Eq. (Q). Details of 
the numerical propagator are given in reference 

In order to go further in the comparison between the 
set of Eqs. (^) and those of STIRAP, we change the repre- 
sentation from the diabatic basis {4>n \ 4>n \ 4>n^} to the 



adiabatic basis of the 3-levels ladder system for each n 
sub-system, {$<_ n) , $jj n) , using the block diagonal 

Rt . . . ® Rn . . . 8 Rn with 



rotation matrix R 



/ sin ip n sin 9 n cos 9 n cos ip n sin 9 n \ 
cos tp n — sin (p n 

sin tp n cos 9 n — sin 9 n cos ip„ cos 9„ 



, (5) 



where the angles are defined by tan# n = flnn/Qrm and 



tan(2^„) = 2 ( ^(ft^) 2 + (tt^) 2 ) /A 



In the new representation, the dynamics of the system 
is followed by the amplitude coefficients {a^\a^\a^}, 
which arc dclocalized and follow the electronic transitions 
between the potentials at different moments of time. The 
new basis is quasi- adiabatic since the transformation ma- 
trix R diagonalizes only each 3x3 sub-system matri- 
ces, while there remain couplings between different sub- 
systems. 

To illustrate the nature of these couplings, we detail 
the analysis for a system of 6-levels obtained by trun- 
cating the expansion of the wave function (Eq. (j|)) to 
only the first two vibrational levels. We consider these 
levels as belonging to two coupled 3-level ladder systems, 
whose Hamiltonian is 



/ fig 1 
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H = 



H21 H 
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n£> 



Q^' -25lu 



(2) 
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(i) 



(G) 



n$ -2(Slj-a) n 22 



(2) 



V n%> 



(2) 
22 



-28u) ) 



and we next change the representation using R = i?i©i? 2 
{R n is given by Eq. (^)). We obtain the quasi-adiabatic 
Hamiltonian 




Ri 1 \ ( H n Vi 2 \(R r 0\ 
%V V ^21 H 22 + ScjTJ\0 R 2 J 



Ri l V l2 R 2 \ 
R^^iRi H {2) +5loI J ' 



(7) 



In this representation each sub-system Hamiltonian 
is diagonal. Neglecting the coupling between sub- 
systems, H® 2 A — R^V^R? ~ 0, and using a counter- 
intuitive sequence of pulses, the Hamiltonian has two 
parallel (independent) dark states as in STIRAP. There- 
fore, the differences between APLIP and STIRAP must 
come from the inter sub-system couplings. Substitut- 
ing cos(6>„) = ilnn/ftnl and sin(6>„) = ftnl/^im where 

Slia = ^(Slml) 2 + (Slml) 2 ^ 2 , we obtain, for our reduced 
6-level Hamiltonian: 
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QA 
12 



cos sin <^2 



(fi 22 n 12 + fi 22 ^12 ) 



"'- P2 coMot 2 ) D^OWl 
„(c) ^'11 J 'l2 "11 "12 > 



cosyi cpsip2 mPlnW i o( 2 'o( 2 H 

T77) S2 12 + S2 11 1Z 12 i 

11 

sinyi sin 1^2 fnWnl 1 ) i o( 2 )o( 2 h 
l"22 "12 "22 "l2 J 



cos /■o( 1 )o( 2 ) 
„(e) v"22 "l2 



sin coWoW _ n( 2 )o( 1 h 
(e) l"22 "l2 "22 "12 ) 



sinyx siny 2 coWoW i o( 2 )o( 2 )\ \ 
<■■>( = ) "-11 12 ~ ^ 11 12 ' 



cos ipi cos 1^2 /-o(l)o(l) o( 2 'o( 2 '"l 
l"22 "12 "22 "l2 J 



n<P2 co(l)o( 2 ) nP'oWl 
~JeT " 11 i,! 12 ~ i2 ll i2 12 J 



cos tp 1 sin go 2 (r>( 1 )o( 1 ) i 
,-j(e) V 11 12 

"11 

sin i^! cos y 2 cot 1 ' n' 1 ) 
TTP, l"22 "l2 



,( 2 ) ( 2 h 



(8) 



Now let us concentrate in the case when both Rabi 
frequencies are equal, [in^i = ^23^2 = ^o- We consider 
the initial state to be the ground vibrational level, |(/>q ) 



which initially correlates with |$g ). Then the only non 
zero couplings between this state and any other quasi- 
adiabatic state are 



= -COS</> 2 (PllP21 



pfxPfx) 



and 



^\H^ ) ) = ^{ P 1 Apfi-pfiPfx) 



n 5i(t)5 2 (t) 



(9) 



(10) 



Equations (^), and (|l0|) separate the molecular contri- 
bution (geometrical factors) from the pulse contribution 
(depending on the pulse amplitude). However it can be 
seen that the APLIP dynamic behavior cannot stem from 
these terms. For instance, in the model of symmetri- 
cally displaced harmonic oscillators the symmetry of the 
system imposes that p\% — P21 and pff = p\±. There- 
fore the quasi-adiabatic initially populated state, I^q 1 "*), 
is not coupled with the rest of the system. So, this state 



is indeed the adiabatic state whose evolution reproduces 
exactly the STIRAP behavior independently from the 
detuning or intensity of the pulses. 

Increasing the number of vibrational levels in the ex- 
pansion (Eq. (^|)) does not introduce any different type 
of couplings and the same conclusions apply. It is very 
simple to generalize the expression for the general 3N- 
level system, in this case, the Hamiltonian matrix terms 
are: 



($<">|jf<M|$( m )) = 

< , 3> (n) I nQ A \(S> ( - m) \ — -On cnitn (r> 12 n 23 - 12 - 23 ^ S 1 (t)S 2 (t) 

($W|ff^|$L m) ) = flosinfatt-pM 



V(piy,(t)) 2 +(pSS2(t)) ! 



For the SDHO model, the geometrical symmetry re- 
quires that: p\ 2 = (^Pl^) = (0i 2) |^i 3) > = Pii, and 
PH = <^Pl4 2) > = (4 2 Vi 3) ) = Pfv Once again, this 
makes any couplings involving the I'&q™' 1 ) states identi- 



cally zero. The system exhibits N parallel (independent) 
trapped states. Therefore, if the initial wave function 
is a coherent superposition of vibrational levels of the 

ground electronic state, \^(x, 0)) = J2n=o a n(0)\(f>iP ( x )), 
a sequence of two coherent counter-intuitive pulses 
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will read this wave function in the adiabatic rep- 
resentation as |*(x,0)> = E^U'inWlSo ^ )). 
that will adiabatically passage at final times to 
\*(x,T)) = E^=o««(0)exp(-i(n-l)^T)|^ 3) (^)) 
(where the phase factor only depends on the eigenval- 
ues of the dressed states |$q™' ) ), i.e. the geometrical or 
Berry phase is zero). At each instant of time only the 
initial levels \<pn^) and the corresponding \<Pn ) levels are 

(3) 

populated, the final population in \(jr n ) being completely 
determined by the initial population, |a n (0)| 2 . The same 
conclusion applies for the final probabilities if the ini- 
tial state is an incoherent sum of vibrational levels. The 
dynamics follows as N independent STIRAP type sys- 
tems and does not reproduce the characteristic features 
of APLIP. 

This result is in flagrant contradiction with the nu- 
merical evidence (direct solution of the TDSE with the 
Hamiltonian of Eq.(l)), in which many of high energy 
vibrational levels of the ground and final potentials arc 
considerably populated (Fig. 1(b)). Although the ana- 



lytic proof derived in this paper is valid for a specific 
"ideal" system, the SDHO model, of course any small 
asymmetry in the model cannot explain the clear dif- 
ference in behavior of STIRAP and APLIP. The differ- 
ence, therefore, must stem from a different source. The 
validity of the quasi-adiabatic Hamiltonian rests upon 
the adiabatic approximation, that is, the neglect of the 

contribution from all the terms coming from R~ 1 R = 
R^Ri^R^^®...® R^Rn, that couple states ®^ 
and <£ > ± l ' > belonging to the same sub-system. In STIRAP, 
the adiabatic condition, Qqt ^> 1 guarantees that these 
terms can be neglected. Surprisingly enough, for very 
large flo and detuning, the APLIP dynamics shows that 
the contribution of these terms cannot be neglected. This 
is in agreement with the numerical results obtained for 
population transfer by shaping LIPs in more general sce- 
narios n Hi ■ 
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